Measurement of the Positive Muon Lifetime and Determination of the Fermi 

Constant to Part-per-Million Precision 
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We report a measurement of the positive muon lifetime to a precision of 1.0 parts per million 
(ppm); it is the most precise particle lifetime ever measured. The experiment used a time-structured, 
low-energy muon beam and a segmented plastic scintillator array to record more than 2 x 10 12 
decays. Two different stopping target configurations were employed in independent data-taking 
periods. The combined results give r M + (MuLan) = 2196980.3(2.2) ps, more than 15 times as precise 
as any previous experiment. The muon lifetime gives the most precise value for the Fermi constant: 
GF(MuLan) = 1.1663788(7) x 10~ 5 GeW 2 (0.6 ppm). It is also used to extract the fi'p singlet 
capture rate, which determines the proton's weak induced pseudoscalar coupling g P . 



A measurement of the positive muon lifetime, t„ 



high precision determines the Fermi constant, 
cording to the relation 
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Here, Aq represents well-known phase space and both 
QED and hadronic radiative corrections and we as- 
sume that Gf is universal for weak interactions. Strictly 
speaking, t /j + determines a muon-decay-specific coupling, 
denoted G M , which could be compared to other Gf de- 
terminations as a test of the standard model 

Prior to 1999, the limitation on the precision of Gf 
was dominated by the uncertainty on Aq. Van Ritber- 
gen and Stuart were the first to compute the second- 
order QED radiative corrections in the massless elec- 
tron limit, reducing the theoretical uncertainty to be- 
low 0.3 ppm and well below the then-current exper- 
imental uncertainty. This development motivated a new 
generation of precision muon lifetime measurements, Mu- 
Lan [|[ and FAST [f| . More recently, Pak and Czarnecki 
extended the result in Q to finite electron mass, which 
shifts the predicted decay rate l/r M by -0.43 ppm; alter- 
natively, it decreases Gf by 0.21 ppm Q. 

In Ref. we reported an 11 ppm measurement of 
t„+ based on a relatively short commissioning run. This 
Letter reports the results from a 100 times larger data 
set, accumulated using the final setup of the experiment. 

The experiment is designed to stop muons in a target 
during a beam-on accumulation interval and measure the 



decay positrons — primarily from the /i + — > e + v e v ll decay 
mode — during a beam-off measurement period. The two 
running periods, R06 and R07, used different targets. 
More than 10 12 decays were recorded in each period. 

The experiment used the 7rE3 beamline at the Paul 
Scherrer Institute (PSI). During the run, positive muons 
from at-rest pion decay near the surface of the production 
target are directed to the experiment through two oppos- 
ing vertical dipole magnets and a series of 15 magnetic 
quadrupole lenses. A velocity-selecting E x B separator 
is tuned to pass muons and reject positrons. A special 
feature of the beamline is a custom, 60-ns switching, 25- 
kV kicker [8|. When energized, the electric field across 
the 120- mm vertical gap and 1500- mm length displaces 
the muon beam by 46 mm at the exit and deflects it by 
45 mrad onto a downstream collimator. The muon flux 
of ~ 10 7 /i/s at the experimental focus is reduced by an 
extinction factor e ~ 10 3 . The kicker is switched by an 
external timing circuit that synchronizes the data collec- 
tion cycles into 5 fj,s kicker-off accumulation periods, Ta, 
followed by 22 [is kicker-on measurement periods, Tm, 
also called "fills." The kicker voltage stability during 
the measuring period, V(t), coupled with the derivative 
de/dV{t) at full voltage, determines the time stability of 
background from this source. 

Muons which arrive at the target are nearly 100% 
polarized. Because of parity violation in the decay, 
positrons are emitted preferentially in the direction of 
the muon spin. Any change in the ensemble-averaged 
residual polarization (P M ) during the measuring period 
can result in an effective time-dependent acceptance if 
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the individual detectors are not symmetric around the 
stopping target and identical in response. For many ma- 
terials, the magnitude of (P^) is reduced during the muon 
stopping process. It can also be reduced by the applica- 
tion of a magnetic field transverse to the spin direction 
and having a field strength sufficient to precess the spin 
rapidly compared to Ta- In practice, a residual longitudi- 
nal polarization (Pl) remains, owing to misalignment of 
the transverse field. The array of detectors is highly and 
symmetrically segmented. The sum of events recorded 
by a detector at angle (8, <fi) with those from an equally 
efficient detector at (180° — 9, 180° + </>) form a decay 
histogram that is immune to changes in (Pp). 

Several changes compared to the setup used in Ref. Q 
involve the muon delivery and the stopping target. A 
200-mm-diameter vacuum beampipe now extends the 
beamline through the detector array (see Fig. [TJ ending 
in a thin mylar window. Muons are stopped in a 21-mm 
horizontal by 10-mm vertical (RMS) spot centered on the 
target disk. The R06 target is a 0.5-mm-thick, 200-mm- 
diameter ferromagnetic AK-3 foil having an internal 
in-plane bulk magnetic field of approximately 0.4 T. The 
R07 target is a 2-mm-thick, 130-mm-diameter crystalline 
quartz (Si02) disk, in which stopped muons form muo- 
nium 90% of the time (lp| . A Halbach arrangement of 
permanent magnets provides a nearly uniform 0.013 T 
field in the plane of the quartz disk. The internal AK- 
3 magnetization or the external magnet field normally 
points left or right in the horizontal direction and the 
muon spin precesses in the vertical plane with periods of 
18 ns and 2.6 ns, respectively, for free muons in AK-3 or 
muonium in quartz. The 10% of muons in diamagnctic 
states in quartz precess with a period of 550 ns, which 
is observed. Both targets can swing open to allow the 
beam to pass unobstructed to a proportional wire cham- 
ber, which monitors the rate and transverse profile at the 
downstream end of the experiment. 

The detector array consists of 170 stacked pairs of 
3-mm-thick plastic scintillators. They are arranged in 
a truncated icosahedron geometry and grouped in 20 
hexagon and 10 pentagon houses. The center of a pair of 
detectors is 383 mm from the target. Each scintillator is 
viewed by either a Photonis (PH) or Electron Tubes (ET) 
29- mm photomultiplier (PMT). On average, 80 photo- 
electrons are registered for each minimum-ionizing parti- 
cle (MIP), which yields a well-defined pulse shape from 
through-going decay positrons. The pulse-height distri- 
bution fit to a Landau function convoluted with a Gaus- 
sian response provides a calibration point for each de- 
tector. The most probable value (MPV) is used as a 
relative calibration of the gain stability versus time in 
the experimental run and, most importantly, versus time 
during the measuring cycle. An LN203 nitrogen laser 
distributes 300-ps-long UV pulses (337 nm) to 24 of the 
detectors and to an externally located reference counter; 
this pulse excites the scintillator similarly to an ionizing 




FIG. 1. Diagram of the experiment. Muons are transported to 
the AK-3 or quartz disk targets through a vacuum beampipe, 
which is lined on the inside with 0.1-mm-thick AK-3 foils 
polarized in the azimuthal direction. 170 pairs of triangular 
scintillator detectors, each read out individually, surround the 
target. The Halbach magnetic ring is used only for the quartz 
disk. A wire chamber (EMC) is placed at the rear to monitor 
the beam profile when the targets are opened. 

particle. The asynchronous pulsing of the laser provides 
a sharply timed common pulse that allows for a relative 
time stability monitor during the measurement period 
and an independent confirmation of gain stability. 

The signals from the 340 PMTs are recorded using 
450 MHz, 8-bit, waveform digitizers (WFDs). Each has 
four independent inputs associated with the inner/outer 
elements of a pair of detectors and the geometrically 
opposite pair. The WFD sampling frequency is set by 
an Agilent E4400 signal generator having stability bet- 
ter than 0.01 ppm/month. The frequency was set to 
±450 ppm of 451 MHz, with the exact value unknown 
to the collaboration and only disclosed after the com- 
pletion of a blinded analysis of the two datasets. The 
"clock tick" (c.t.) time units have the conversion 1 c.t. 

2.2 ns. Normally 24 ADC samples (53 ns) comprise 
a full waveform; they are recorded when an input signal 
exceeds a set threshold, with four to eight samples pre- 
ceding the trigger point. The digitization is extended if 
the threshold is exceeded at the 24th sample. For the 
average pulse, the full width at 20% maximum is 9 ns. 

The 130 TB of raw waveforms are first converted into 
lists of valid "hits" having characteristic times and ener- 
gies for each detector. The waveforms are fit to pulse- 
shape templates, prepared from a set of low-rate events; 
the time of a decay is defined as the peak of the pulse 
shape. In more than 99.9% of the cases, a single pulse 
exists on a waveform and it can be identified reliably. 
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A distribution of pulse energies is made for each detec- 
tor and a run-by-run dependent pulse threshold is estab- 
lished at the minimum in the distribution between the 
low-amplitude PMT noise and the MIP peak (45 - 55% 
of the MIP peak amplitude) . Rare multiple-pulse wave- 
forms are more difficult to characterize, requiring an it- 
erative approach that adds additional pulses as needed 
and/or introduces more sophisticated fitting methods. 
When two pulses are separated by more than 3 c.t., 
the procedure works well. For shorter separation times, 
only one hit is reconstructed and its time is set to the 
energy-weighted average of the ADC samples. For re- 
solved pulses, an "artificial deadtime" (ADT) can still be 
applied on a per-detector basis, eliminating pulses when 
a minimum time separates sequential hits in the same 
scintillator. The ADT used in the analysis varies from 5 
to 68 c.t. and represents an important diagnostic of the 
pileup correction procedure. 

Histograms are filled with inner-outer coincident 
events for each detector pair; the coincidence window 
width is set equal to the ADT. The muon lifetime is ob- 
tained from a fit to the sum of the 170 individual his- 
tograms using the function 

F(t)=Aexp(-t/Tp) + B, (2) 

where B accounts for the flat background. Under 
our running conditions, uncorrected pileup would shift 
the fitted lifetime from its true value by 5t,,/ ADT = 
24 ppm/c.t. Pileup is corrected using a statistical pro- 
cedure based on the data itself. When a hit is observed 
at time U, in fill j, an interval between U and tj+ ADT 
is searched in fill j + 1. If a hit is observed in this inter- 
val, its time is recorded in a separate histogram, which 
is then added back into the original decay histogram. 
The process is repeated for higher-order pileup, and it in- 
cludes fine adjustments to account for time jitter between 
inner-outer coincidences, and for uncorrelated single hits 
and other distortions. A Monte-Carlo study verified the 
pileup correction procedure with input hits occurring at 
actual coincidence, singles, and background rates, and 
with the measured time jitter in inner-outer hits. It was 
also tested at much higher rates where the correction is 
correspondingly larger. Pileup-related lifetime shifts are 
reduced to Jr M /ADT = +0.018 ppm/c.t. A linear extrap- 
olation to zero ADT gives our quoted value. The overall 
pileup systematic uncertainty is estimated to be less than 
0.2 ppm. An alternative approach is to pileup correction 
is to add an explicit lst-order pileup term, exp(— 2t/r iU ) ) 
to the fit function in Eq. [5J This yields the same lifetime 
as our correction procedure, but at the price of doubling 
the uncertainty on r^, owing to the correlation of the 
additional term in the fit. 

Because the residual polarization in the AK-3 target is 
small enough that it is not distinguishable in individual 
time spectra, Eq. [5] describes the R06 data well. A fit 
to the pileup-corrected event histogram, summed over 



all detector pairs, gives a x 2 /dof = 1.03 ± 0.04 for a fit 
start time 220 ns into the measurement period and an 
ADT of 6 c.t. Although no <f> dependence is observed for 
individual fits, a 12 ppm slope exists in the distribution 
of lifetime versus cos(0), which is canceled to better than 
0.10 ppm by the symmetry of the detector array and the 
accuracy of the alignment. The slope is consistent with 
the relaxation of a spin component along the beam axis. 
The lifetimes of the 85 opposite detector pairs, which are 
first summed, then fit, is flat, when plotted from 9 = 
to 90 degrees, and the distribution of results is Gaussian 
having a normalized width cr/mean = 0.90 ± 0.08, which 
is in acceptable agreement with 1.0. 

In contrast, the R07 quartz data are complicated by 
the relatively slow precession and relaxation of the dia- 
magnetic muons. The amplitude of the precession sig- 
nal for individual detectors can be as large as 0.1% and 
the longitudinal polarization can be as large as 0.15%. 
A different analysis procedure, which is insensitive to 
small and unknown efficiency differences, is used. Each 
detector-pair histogram is first fit using Eq.[2]with A mul- 
tiplied by [1 + P 2 -exp (— t/r 2 ) sm(uit + </>)]; here Pi and r 2 
are the transverse polarization and relaxation time con- 
stant, respectively. These fits result in 170 "effective" 
muon disappearance rates, which differ from the true 
muon decay rate because of the relaxation of the longi- 
tudinal spin component. Next, the 170 effective lifetimes 
are fit to extract the true muon lifetime using a func- 
tion that accounts for the position of detectors relative 
to the longitudinal polarization direction. This fit takes 
advantage of the smooth geometrical-based evolution of 
the relaxation effect. The leading systematic uncertainty 
of 0.20 ppm for this procedure is due to the uncertainty 
of the beam position on the target. The lifetime obtained 
is in agreement (to 0.3 ppm) with the result from a fit to 
the simple sum of detector-pair histograms using Eq. [2l 

Instability of the combined detector and electronics re- 
sponse, or "gain," during the measurement period is a 
possible source for a systematic error. It will result in 
gaining or losing hits, which can shift the fitted lifetime. 
Four effects were identified: (i) When triggered to begin 
the recording of fills, a fanout unit induced a 1.3 x 10~ 4 
gain oscillation in each WFD. The effect, which is repro- 
ducible in bench tests, disappears in ~ 4 fis. An empir- 
ical function is used to characterize and correct for the 
oscillation prior to the final fits, (ii) A short-time-scale 
gain shift occurs when one hit closely follows a preced- 
ing hit. The second hit is affected by the remnant long- 
decay-time scintillator light from the first pulse and by 
the intrinsic recovery time of the PMT and voltage di- 
vider. A ~ 1% increase in detector gain exists following 
a pulse, which persists for 50 or 500 ns in PH- or ET- 
instrumented detectors, respectively, (iii) A long-time- 
scale change also exists that differs by PMT type. The 
PH-tube gains increase by 1.8 x 10 -4 and the ET-tube 
gains decrease by 5 x 10~ 4 , over the 22 /is measurement 
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period, (iv) When below-threshold pulses fall either in 
the peak or pedestal regions of a trigger pulse, they can 
either raise or lower the trigger pulse amplitude, respec- 
tively. Strictly speaking, this is a pileup effect, which is 
not accounted for in the pileup-correction procedure. 

These four phenomena are interrelated and give rise 
to an overall variation of 3 x 10~ 4 in the gain during 
the fill. This variation was separately evaluated for PH- 
and ET-instrumented detectors using the MPV versus 
time-in-fill, and was used, after normalization and prior 
to the final lifetime fits, to correct for the resulting time 
spectrum distortion. After all corrections are applied, 
the shift in lifetime compared to uncorrected spectra is 
+0.50 ± 0.25 ppm. The procedure was tested by accu- 
mulating a high-threshold spectrum, which amplifies the 
time spectrum distortion by roughly a factor of 60. The 
lifetimes derived from normal- and high-threshold set- 
tings, following the correction procedures, are in good 
agreement. 

The complete summary of systematic uncertainties is 
given in Table HI including small effects related to muon 
stops upstream of the target, the short-term and long- 
term time response stability based on the laser system, 
and the uncertainty on the absolute clock frequency. 



TABLE I. Systematic and statistical uncertainties in ppm. 
The errors in different rows of the table are not correlated to 
each other. If only one error appears in a given row, the effect 
is 100% correlated between the two run periods. 



Effect uncertainty in ppm 


R06 R07 


Kicker stability 


0.20 0.07 


Spin precession / relaxation 


0.10 0.20 


Pileup 


0.20 


Gain stability 


0.25 


Upstream muon stops 


0.10 


Timing stability 


0.12 


Clock calibration 


0.03 


Total systematic 


0.42 0.42 


Statistical uncertainty 


1.14 1.68 



The stability of the result versus start time of the fit is 
a powerful collective diagnostic because pileup, gain sta- 
bility and spin effects all might exhibit time dependence. 
For both R06 and R07, the lifetime does not depend on 
the fit start time, apart from the statistically allowed 
variation. Further, it does not depend on run number or 
magnetic field orientation. 

The final results for the two running periods are in 
excellent agreement: 



t m (R06) = 2196979.9 ± 2.5 ± 0.9 ps, 
t m (R07) = 2196981.2 ±3.7 ±0.9 ps. 



(3) 




Balandin - 1974 
Giovanetti - 1984 
Bardin - 1984 
Chitwood - 2007 
Barczyk - 2008 
MuLan - R06 
MuLan - R07 



FIG. 2. Lifetime measurement summary. The MuLan R06 
and R07 results are plotted separately to illustrate the con- 
sistency. The vertical shaded band is centered on the MuLan 
weighted average with a width equal to the combined uncer- 
tainty. 



at the ppm level, the expectation that the lifetime of 
bound muonium does not differ appreciably from the free 
lifetime [ll|. Combined we obtain 

r M (MuLan) = 2196980.3 ± 2.2 ps (1.0 ppm), (4) 

which is in agreement with our previous measurement [J] . 
The error is the quadrature average of statistical and sys- 
tematic errors, where the full error matrix calculation, in- 
cluding all correlations, is used to combine uncertainties. 
The MuLan result is more than 15 times as precise as 
any other individual measurement and consequently 
dominates the world avera ge. Our result lies 2.5 a below 
the current PDG average 12j . Figure [H shows the re- 
cent history of measurements, together with the MuLan 
average. 

Our value for 
tion of the Fermi constant: 



, ,,-r leads to the most precise determina- 



Gf (MuLan) = 1.1663788(7) x 10 _t> GeV^ (0.6 ppm). 

(5) 

The positive muon lifetime is also used to obtain ordinary 
muon capture rates in hydrogen [l3| or deuterium [l4| 



by the lifetime difference method, T c 



1/V - l/r„ 



These capture rates determine hadronic quantities as dis- 
cussed in |15[. For example, our new result lowers the 



fi~p capture rate used in Ref. |13[ by 8 s , and thus 
shifts g P upward to even better agreement with theory. 
Finally, the improved precision reduces the uncer- 
tainty in the determination of muon capture in [l3j and 
future efforts to below 0.5 s -1 . 

We thank the PSI staff, especially D. Renkcr, K. Deit- 
ers, and M. Hildebrandt; M. Barnes and G. Wait from 
TRIUMF for the kicker design, the NCSA for enabling 
and supporting the data analysis effort, and the U.S. Na- 
tional Science Foundation for their financial support. 



Here, the first errors are statistical and the second sys- 
tematic. The comparison between R06 and R07 affirms, 
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